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I propose a variation of the standard worm-like chain model to account for internal order param- 
eter (helix/coil) fields on the polymer chain. This internal order parameter field infiuences polymer 
conformational statistics by locally modifying the persistence length of the chain. Using this model, 
I make predictions for the bending and stretching response of an alpha-heliacal domain of a protein. 
In particular, I show that alpha-helical protein domains will buckle under applied torque. This 
highly nonlinear elastic behavior may be important in the understanding of allosteric control of 
biochemical pathways. 



PACS numbers: 87.10.-|-e,87.14.Ee,82.35.Lr 

The principal coarse-grained model of a semiflexible 
polymer is the worm-like chain (WLC). It describes the 
single-chain polymer statistics in terms of a quadratic 
Hamiltonian that associates an energy cost with chain 
curvature by introducing a bending modulus k. In ther- 
mal equilibrium the effect of this bending energy is to 
enhance the statistical weight of straight polymer confor- 
mations on a length scale k (/cbT = 1) called the thermal 
persistence lengthTj. The effect of this bending energy 
at longer length scales is trivial and the single-chain sta- 
tistical mechanics at those longer scales can once again 
be described by standard Gaussian statistics (ignoring 
for the moment the effect of solvent quality) with an ad- 
justed Kuhn length. 

In this letter I propose a coarse-grained polymer model 
along the lines of the WLC designed to describe polymers 
with internal degrees of freedom. In particular, I seek to 
describe the statistical mechanics of alpha-helix forming 
domains of proteins(^] in an analytically tractable way. 
By using such a coarse-grained model for the (nonlinear) 
elastic properties of a single protein domain, one can at- 
tempt to construct a low-dimensional description of pro- 
tein mechanics (i.e. conformational change) in terms of 
a small set of locally interacting subdomains of definite 
secondary structure. In particular, a generic consequence 
of these simple models is that due to the highly nonlin- 
ear nature of the effective elasticity of an alpha helix, 
these structures are bistable under an externally applied 
torque. Their mechanical bistability should naturally de- 
scribe the allosteric control of enzymes 0- 

The internal degrees of freedom are two-state Ising 
variables that represent the presence or absence of hy- 
drogen bonding between the helix turns of the molecule 
in its native secondary structure. The presence of these 
internal state variables - hydrogen bonding - couples to 
the conformational degrees of freedom of the molecule 
through the thermal persistence length. In the hydrogen- 
bonded state the polymer is stiffer having a larger bend- 
ing modulus since bending the molecule involves the 
stretching and compression of the hydrogen bonds. 

The internal state variable describing the pres- 



ence/absence of native secondary structure and conse- 
quently local hydrogen bonding is also controlled by its 
own energetics in addition to the aforementioned cou- 
pling of these internal degrees of freedom to the con- 
formational degrees of freedom of the space-curve repre- 
senting the polymer chain. The Hamiltonian for these in- 
ternal degrees of freedom is a one-dimensional, ferromag- 
netic Ising model that in the current context has been re- 
ferred to as the "helix-coil" (HC) model Q ■ This model 
has been previous employed to study a variety of inter- 
nal state transformations of biological macromolecules 
such as the unzipping of double stranded DNA, and the 
opening-up of alpha-heliacal proteins under applied ten- 
sion [3. However, this model has never been combined 
with the WLC model of the conformational degrees of 
freedom of semi-flexible polymers in order to explore 
in detail the coupling of these internal degrees of free- 
dom along the protein backbone with both bending and 
stretching degrees of freedom. The new model, the helix- 
coil worm-like chain (HCWLC) combines the internal 
of helix-coil degrees of freedom with the conformational 
states of the semiflexible chain. 

The HCWLC Hamiltonian is a function of discrete set 
unit tangent vectors ti, i = l,...iV-fl and Ising vari- 
ables Si = ±l,i = l,...iV + l. In this coarse-grained 
description of the polymer, each segment of the chain has 
two variables associated with it - an Ising variable de- 
scribing whether that segment of the chain maintains its 
native (alpha-hehcal) secondary structure (s — +1) or 
not (s = —1) and a local tangent vector t lying along the 
averaged direction of that segment of the chain. See fig- 
ure one. Thus, I do not consider more than two internal 
states of the polymer, but clearly such a generalization 
(creating a q-state Potts model along the Id chain) is 
possible. In addition, I do not consider here a torsional 
degree of freedom along the chain. In order to reduce 
the calculation to its most fundamental level, I will also 
analyze the model in two dimensions so that the tan- 
gent vectors are equivalent to an angular variable. Three 
dimensional calculations follow similarly [||. Finally, it 
is important to note that segments of the chain so de- 
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scribed represent more than a single monomer since these 
segments must he assigned a particular, local secondary 
structure. Based on these considerations, the Hamilto- 
nian is given by 

N , N+1 

1=0 2 = 

N 

+ ^k(s,)[1-cos(0,+i -0,)], (1) 

1=0 

where the bending modulus of the chain, k(s) takes the 
form: 



k{s) 



if s = +1 
K<g if s = — 1 



(2) 



There are four parameters with the dimensions of en- 
ergy in Eqs. I1I2I The first appearing in the Hamiltonian 
is ew that represents the energy cost of a domain wall 
in the internal Ising-like variable. Physically this energy 
corresponds to the free energy cost associated with ei- 
ther starting or ending a sequence of alpha-helical seg- 
ments of the polymer and thus reflects a combination 
of hydrogen-bonding between consecutive turns of the 
alpha-helix and the loss of local entropic configurational 
free energy due to the adoption of helical structure. In 
the standard helix-coil terminology, this energy is the 
logarithm of the "chain cooperativity" parameter. The 
second parameter in the model is h, which appears as a 
fictitious magnetic field favoring the spin up {i.e. alpha- 
helical secondary structure) configuration of the chain. 
In the current interpretation of the HCWLC model, this 
term represents effect of the local, solvent-polymer in- 
teractions which provide a thermodynamic driving force 
toward its native secondary structure (alpha-helix) per 
unit length. The final two energies enter in the last term 
of Eq. n as seen in Eq. [21 the local bending modulus of 
the chain takes one of two values depending on whether 
that segment of the chain has its native structure (s = 1) 
or not (s = ~1)- Fundamental to the basic phenomenol- 
ogy discussed below, the bending energy of the chain is 
higher in the alpha helical state due to the presence of 
hydrogen bonding so that generically ^ k^. 

One can now compute the response of the chain to ap- 
plied torques and forces in equilibrium - quantities can 
be directly measured in single-molecule manipulation ex- 
periments. I consider the response of molecule to applied 
torques which constrain the relative angle between the 
two end tangent vectors of the chain to be a. Given that 
constraint, the partition fmiction of the polymer is given 
by 

N+1 „ 

2(a) =n E de.a,^-"'^ 5{Oo)5{ON+^-a), (3) 

i=0 Si=±l-' 

where the sum is over all angles and internal states of the 
chain while the two delta functions fix the initial tangent 



FIG. 1: The upper figure represents the description of a poly- 
mer with secondary structure. The open circles represent the 
internal state variable Si characterizing secondary structure 
while vectors are the local tangents to the space curve of the 
polymer. The lower figure shows that space curve (consistent 
with boundary conditions introduced in Eq. I^J but does not 
explicitly show the secondary structure. 



vector to lie along the x axis (removing the trivial, degree 
freedom associated with uniform rotations of the system) 
and the final tangent vector to make an angle a with the 
initial one thereby enforcing a particular total bend in 
the polymer chain. 

The main qualitative result of this calculation is that, 
due to the coupling of the conformation and internal de- 
grees of freedom of the chain, the polymer is bistable 
under applied torque. Upon increasing the angle a, one 
finds that the necessary torque increases linearly for small 
angles. Upon reaching a critical angle, however, this 
torque drops dramatically indicating the nucleation of 
a soft (i.e. short persistence length) region of the chain. 
This new instability, which is due to the coupling of the 
chain persistence length to its state of secondary struc- 
ture, should have a clear experimental signature in single 
protein manipulation experiments. 

I evaluate the partition function by a transfer matrix 
technique^. Using the identity that relates the expo- 
nentiated cosine to a sum of modified Bessel functions of 
integer order, one can recast the partition sum in terms of 
a set of variables conjugate to the angles 9i, i = 1, . . . N 
(quantized angular momenta rm of a particle on the unit 
circle)that also form a complete set of states. Using this 
completeness the partition sum in Eq. |31 can be written 
as 

2(") = X! Vsn+i{so,(^q\T^^^\sn+i,On+i ^ a) 

SO,SiV + l = ±l 

(4) 

where initial ((|) and final (|)) states of the chain are 
two-component objects representing both the secondary 
structure of the these segments (s) and the angle of their 
tangent vector, which has been determined by the delta 
functions appearing in Eq. |3| There is a one-to-one map- 
ping between partition function of a chain of A'^ -I- 1 seg- 
ments having its initial and final tangents at angles and 
a respectively and the quantum transition amplitude of 
a particle moving on the unit circle from angle to angle 
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a in N time slices. Local secondary structure acts like an 
internal two-state variable analogous to the spin state of 
a spin— 1/2 particle. The transfer matrix (or time evolu- 
tion operator) T is diagonal in a momentum taking the 
form of a 2 X 2 matrix acting on the "spinor" space of 
secondary structure so that: 



we write the partition function as 



{m, s \T\m' , s') — 

/m(K»)e-''» /m(K»)e~'^» 
-h 



(5) 



The term rjs ~ 5s,-ie^^ in Eq. 01 is a book-keeping de- 
vice required to include the free energy cost the N + 1 
segment of the chain not being in its native state. This 
asymmetry between the two ends of the chain is an arte- 
fact of the model that is due to the fact that the bending 
modulus associated with the v}^ and [n -\- 1)*"^ tangent 
vectors depends on the spin s„ only. Thus the quantum 
analogue of this problem described by the above Haniilto- 
nian is non-Hermitian and consequently asymmetric with 
respect to time reversal. 

The transfer matrix defined in Eq. [S] can now be triv- 
ially diagonalized in the spinor space of secondary struc- 
ture by a similarity transformation. The eigenvalues of 
the transfer matrix are given by 



Ai,2(m) = ^^^^^ h+z^±J{l- z^f + a^z. 



(6) 

where = exp(— is the scalar operator 

that plays the role of the transfer matrix for the all al- 
pha helical, WLC, while z^ — exp(— /i)[x)„i(k^)/cj,„(k3>) 
is the ratio of the fugacity of a random coil (non-native 
state) segment to the alpha-helical fugacity both evalu- 
ated at angular momentum m. Finally, (3 is the expo- 
nentiated free-energy cost of introducing a domain wall 
in the spinor variables representing a change in the local 
secondary structure: /3 = exp(2 — e^), where the first 
term reflects the increase in entropy associated with the 
appearance of one domain wall, and the second term is 
the energetic cost of such a structure. 

The diagonalization of the transfer matrix has been 
effected by the similarity transformation: 7"^+^ — > 
. T^^^U. In the thermodynamic limit N — > cx), 
one may impose cyclic boundary conditions and then use 
the properties of the trace to eliminate the matrix U with 
its inverse. Here, it is important to keep finite to study 
the effect of the degree of polymerization N on the phe- 
nomenology of the model. In addition, alpha-helical seg- 
ments of proteins typically contain no more than order 
ten alpha-helical turns so the quantitative application of 
the model to real protein mechanics demands an accu- 
rate description of the system at these small degrees of 
polymerization. Using the eigenvalues defined in Eq. 



Z(a) 



oo , 
^ g^am I _ 

m— — OO ^ 

Af (m) ~ Af (m) 
2^(l-z„)2+/52^„ 



Af(m) + A^(m)] 



-/?])}. (7) 



One can gain insight into this expression in the limit of 
high chain cooperativity where ^ 1, X2{'m) <C 1 by 
expanding in powers of p. By collecting terms in this 
way Z = Z*^"^ -I- /JZ^^^ we reorganize the partition 

sum into a term counting all conformational states of the 
all alpha- helix and all random coil WLC (Z^"^) plus all 
states of the polymer with one domain wall (Z'^^) etc. 

The second term of Eq. [7| is due to the nonperiodic 
boundary conditions and is the sole source of such odd- 
order P terms. Viewing the appearance of random coil 
segments of the chain as a nucleation problem, one may 
refer to terms as representing a form of heterogeneous 
nucleation; in this limit the domain wall energy is so high 
that "melting" the alpha-helical structure from the chain 
ends is energetically favorable. For longer chain lengths 
and/or for smaller domain wall energies, melting inter- 
nal segments becomes favorable due to increased transla- 
tional entropy associated with the placement of random- 
coil segment compared to the energetic cost of a second 
domain wall. This transition to internal, or homogeneous 
nucleation of random-coil segments should occur at chain 
lengths where logiV ~ e^. For alpha- helical domains of 
e.g. six helical turns, this transition between the two 
forms of melting occurs at domain wall energies around 
1.5fcBT. 

Based on the partition function given in Eq. |7| we can 
compute all statistical properties of the chain for which 
the directions of the ends are constrained. We two most 
interesting quantities are: (1) the mean torque required 
to fix a given angular bend of the chain T(a) and (2) the 
fraction of the chain in the non- native (random coil) state 
as function of the same angle M(a). These quantities are 
given respectively by 



T{a) 
M{a) 



dlnZ 
da 

1 91nZ 
~N dh ■ 



(8) 
(9) 



Both of these quantities can be expressed as sums over 
one remaining angular momentum by taking the appro- 
priate derivatives of Eq. \7\ These infinite sums are exact 
results. The torque angle curve (approximated by a par- 
tial sum) is shown in figure |21 Keeping the first fifteen 
terms gives this curve with an accuracy of one part in 
10^. 

The principal qualitative feature of these torque angle 
curves is that for small angles the a-helical chain bends 
as a simple, semiflexible object; at a critical angle, how- 
ever, the free energy of the chain is minimized by the 
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FIG. 2: The upper (dashed) curve shows the torque r (mea- 
sured in units of ksT vs. angle a curve for the model with 
13 = 10"*. The lower (solid) curve is the equivalent plot for 
l3 = 10"^. For both curves, N = 15, = 100, k« = 1, 
ft = 3. 

creation of a random-coil segment. The free energy in- 
crease associated with the introduction of the thermo- 
dynamically unfavored random-coil segment is offset by 
the decrease in curvature energy due to the fact that the 
chain can now bend at the more flexible link thereby cre- 
ated. It may be noted that for the physically relevant 
quantities used in figure |21 the critical torque is on the 
scale 5 — lOfceT « 20 — 40pN ■ nm; such torques are 
clearly achievable through e.g. molecular motors acting 
over length scales of merely a few nanometers. One can 
observe the creation of random-coil segments by com- 
puting the magnetization, or fraction of the chain with 
native secondary structure. At the buckling angle (the 
maxima of the curves in figure [2Jl, M ~ 0{1/N) demon- 
strating the nucleation a single weak link. The complete 
set of curves will be shown later j^]. 

Another testable experimental prediction involves the 
force vs. extension curves of such a molecule. Such cal- 
culations have been performed for the WLC model0,l3- 
While the torque angle curves can be simply computed in 
an exact manner due to the fact that the angle variables 
form a natural representation of the semiflexible chain, 
the force extension behavior of the model is somewhat 
less transparent. These calculations can, however, be 
preformed perturbatively in the strength of the applied 
force. The simple structure of the perturbation series 
(shown below) allows the efficient computation of terms 
to arbitrary order in the applied force. To introduce ap- 



plied forces, we include a potential in the Hamiltonian H 
from Eq. ^o that 

N+l 

if — >i/-/^7(s.)cos0, (10) 

i=0 

where 7(5) is the length of a segment, which depends 
on the local secondary structure: 7(+l) < 7(— 1) and 
/ is the applied force stretching the end of the chain 
in the x direction. The exact evaluation of the partition 
sum is now quite complex since the complete set of states 
needed to insert between time-evolution operators is now 
the wave function of the quantum pendulum ( Matheau 
functions). The extensional response of the chain to an 
applied force, however, can be computed perturbatively 
in powers of /: 

N+l 

Z{aJ) = Z{aJ ^ 0) + fJ2{lis^) cos 9,)o+ (11) 

i=0 

^ N+l 

+ 2/^ (7(s0cos6',7(sj)cos6lj)o + --- 

where the thermal averages (■)o are taken with respect 
to the Hamiltonian with f = The computation 
of these averages is straightforward using the transfer 
matrix method. 

HCWLC theory is the simplest extension of the WLC 
admitting a coupling of the chain conformational degrees 
of freedom to internal state variables. It replaces one per- 
sistence length of the WLC with four energy scales that 
must be determined via experiment, however, it gener- 
ically creates bistable configurations of the chain under 
applied torque. The bistability of the individual subdo- 
mains of proteins provides a framework for understanding 
protein conformational change. 
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